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The thermodynamic properties of proton rich systems are explored in a mean field approach which 
is generated from a Skyrme interaction. The addition of Coulomb interactions result in asymmetries 
which modify the chemical and mechanical instability of the system and its equilibrium properties. 
These properties are studied for systems with proton fraction y on the proton richer side of the 
valley of /3-stability as well as the neutron rich side. Coulomb induced instabilities lead to proton 
C**) ■ diffusion processes on the proton richer side and also large asymmetries in chemical and mechanical 

instabilities and coexistence curves. Considering the whole range of < y < 1, we can study how 
the symmetry about y = 1/2 is broken by asymmetric interaction and we can also explicitly show 
that the role between proton and neutron is exchanged around yE- It is shown that there are two 
(3JT)' asymmetric coexistence surfaces in (y,P,T) space, one for y < yE and another for y > %)e and 

^ t touching each other at yE- These asymmetries in instabilities show up as new branches, one for 

y < yE and one for y > yE, and thus form a closed loop in pressure versus p for both chemical 
instability and coexistence regions. The branch of y > 1/2 > yE was not previously investigated 
since only the y < 1/2 region is usually considered. In our simplified model, mechanical instability 
is still symmetric around a point yE 1/2 even with Coulomb forces present. 
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Understanding properties of nuclei with large asymmetry in the proton neutron ratio has had a renewed interest 
for several reasons. For example, future experiments done with rare isotope accelerators (RIA) can explore properties 
of nuclei far from the valley of /3-stability on both the neutron and proton rich sides of it. These studies would be 
useful for understanding the production of the elements and neutron stars in the area of nuclear astrophysics. Current 
interest stems from the liquid-gas phase transition which is experimentally studied using medium energy heavy 
ion reactions. A primary goal of heavy ion reactions is to explore the phase diagram of hadronic matter both at low 
1 and high density and temperature looking for signals of phase transitions. Another reason for studying asymmetric 
two component nuclear systems is that they offer a unique example of the thermodynamic and statistical properties 
of binary systems. Here, the nuclear phase diagram is governed by nuclear volume, surface, symmetry energy and 
asymmetric Coulomb effects. The interplay of these on the phase structure has been of interest in past studies 
Examples of other two component systems are binary alloys and liquid 3 He. For 3 He the two components are spin up 
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and spin down fluids and important new phase structure appears within the superfluid phase itself below 3 mK. 

In this paper we extended our previous study 0iQ of phase structure of nuclear systems below the liquid-gas critical 
point to include the entire range of the proton fraction y — Z/{N + Z) from to 1. Our previous study focussed 
on < y < 1/2. For systems with no asymmetric Coulomb effects, a symmetry exists about the y = 1/2 point. 
This symmetry is broken by the asymmetric Coulomb term. Due to this symmetry, usually only y < 1/2 side was 
considered and thus such studies would miss various symmetry breaking features and the existence of another pair of 
coexistence points at large y. Some preliminary initial new findings were presented in Ref.@,ll] in a restricted region 
of y < 1/2 which are further elaborated on here and also extended to the whole range of < y < 1. In particular we 
will focus on the proton rich side of the valley of /3-stability and study the chemical and mechanical instabilities and 
equilibrium properties of these systems. Several previous mean field studies of two component systems were carried 
out without the Coulomb force. The pioneering work of Miiller and Serot focussed on the importance of the nuclear 
symmetry energy on the equilibrium structure of liquid-gas phase transition. V. Baran etal [j| investigated instability 
properties of such systems and pointed some new properties of the liquid-gas phase transition associated with chemical 
instability and isospin distillation based on symmetry energy considerations. Calculations were extended by Colonna 
etal to finite nuclei including surface and Coulomb effects. Their calculation focussed on the instability modes and 
their results showed that an octupolar mode dominates. Pawlowski [j| has also recently considered the importance 
of Coulomb and surface effects on the equilibrium phase structure of the liquid-gas phase transition and found them 
to be important. An extensive summary of isospin effects in heavy-ion collisions at medium energy can be found in 
the book of Bao-An Li and Schroder Q . 

In our study we will use a simplified Skyrme interaction to generate the nuclear long range attraction, short range 
repulsion and also the nuclear symmetry energy. We include a simplified Coulomb energy which generates a proton- 
neutron asymmetry. The resulting interaction energy density is 

U(p) = -a oP 2 + a 3 p 3 + a s (2y - if p 2 + Cy 2 p 2 (1) 
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The do term is the long range attraction which is quadratic in the density p. Shorter range repulsion is the 03 term and 
we take its density dependence to be p 3 . The nuclear symmetry energy is the a${2y — l) 2 term with as = |(| +Xo)aQ 
and has a p 2 dependence. The asymmetric Coulomb term is Cy 2 p 2 . The C, which in principle depends on nuclear 
size R, is taken to be constant. For this choice, we can consider R to be an effective range of the Coulomb interaction. 
Then, in this simplified limit, the Coulomb interaction depends on the square of the proton density. Our choice for the 
Coulomb energy and nuclear interaction energy greately simplifies the full complexity of the analysis and results in 
analytical solutions for instability boundaries. By using a somewhat simplified imteraction, we hope to qualitatively 
illustrate the main physical phenomena associated with phase transitions in two component systems which include 
both symmetry and asymmetric terms. We also believe that the main results will not change qualitatively, nor in any 
major quantitative way when more complicated dependences are included. 

A kinetic energy density is added to U(p). Since we will present results for moderate density and temperatures of 
the order of 10 MeV, we will work to first order in the degeneracy corrections. In this limit the pressure is 

p = > T - {" - 57ll) ' 2 + 2 '-" 3 + (§<5 + + w£) (2 » - 1)V + c » v (2) 

The a = -§i and a 3 = ±t 3 with t = -1089 MeV fm 3 and t 3 = 17480.4 MeV fm 6 as shown in Ref. 2]. The 7 = 2 
is the spin degeneracy factor and A = hj ' (^irmT) 1 / 2 . The proton and neutron chemical potentials are 

p q = T\n[(\ 3 / 1 )(l±(2y-l))p/2}-2a oP + 3a 3 p 2 

±(4/3)(l/2 + x Q )a Q (2y - l)p + (1 ± \)Cyp + (T/2\/2)(A 3 / 7 )[l ± (2y - l)]p/2 (3) 

with the upper sign for protons and the lower sign for neutrons. These three equations determine the equilibrium 
phase structure and the mechanical and chemical instability regions of our simplified two component nuclear system. 
The phase structure is a three dimensional surface in y, T and P. At y = 1/2, symmetric systems with no Coulomb 
interaction, one has the familiar P, T curve of Maxwell pressure versus T for equilibrium state which ends at the 
critical point. For fixed T, the equilibrium surface intersects the fixed T plane to form loops in P versus y. These 
loops represent the proton fraction in the liquid and gas phases which arc different except at a point called yE which 
is not 1/2 for the case with Coulomb interaction. These equilibrium loops end at the critical points (y c (T), P C (T)) at 
fixed T or (P c (y), T c (y)) at fixed y having the same densities both in liquid and gas phases. The point with lowest y 
value on this loop is the maximally asymmetric point at a fixed T. Another maximal asymmetric point is at y > yE 
with largest y value. The equilibrium surface is generated from the chemical potentials p q {y,T 1 P) by requiring at 
fixed P and T the neutron chemical potential in liquid and gas phase to be equal and proton chemical potential in 
these two phases to be equal at the same two values of y. Geometrically, this corresponds to the construction of a 
rectangular box in a plot of p, versus y at fixed P and T where the horizontal lines are the same neutron chemical 
potentials and the same proton chemical potentials and the vertical lines are the two values of y where each chemical 
potential in liquid and gas phases are equal. The gas phase will be more asymmetric and have lower y for y < yE and 
larger y for y > yE than the liquid phase. The liquid phase is more symmetric because of the symmetry energy term 
and has the value of y closer to yE- Figures ^ an d [21 shows various properties associated with the equilibrium surfaces 
and behavior of the chemical potentials. 

The point where the proton fraction is equal in both gas and liquid phase can be obtained from the condition 
dP I dy = and this point is given by 

An important feature of the Coulomb force is to shift yE away from its value of 1/2 when C = to a value closer 
to the valley of /^stability found from d(E / A) / dp = P = at zero T. This shift is important for understanding the 
phenomena of isospin fractionization which favors more symmetric liquids and less symmetric gases. 

FigD shows the effect of various terms on the coexistence curve obtained from the intersection of the coexistence 
surface with a fixed T = 10 MeV plane. Here we have used the same parameters as used in Ref. 0,0 except for using 
a more realistic smaller symmetry energy (xo = —1/6). The curve has two symmetric loops which are connected with 
minimum P at y — 1/2 when Coulomb forces are turned off and the loops are larger for a large symmetry energy 
(dotted line with xo = 1/2) compared to smaller symmetry energy (dashed line with xq = —1/6). A comparison 
of these two curves show that the symmetry energy has essentially no effect on the minimum point at y = 1/2 
(symmetry energy is zero at this y) while the maximum points (which are not at y = 1/2) has larger pressure for a 
larger symmetry energy. The dash-dotted line with minimum at y — 1/2 includes a surface energy or surface tension 
term besides a symmetry energy term whose value is the same as that for the dashed line. By comparing these two 
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FIG. 1: Pressure P versus proton fraction y for coexistence loop at T = 10 MeV. The solid line is for the case with Coulomb and 
surface effects. The dash-dotted line has surface, but no Coulomb terms. Coulomb effects are shown in the dash-dot-dot-dotted 
line with no surface terms. The dashed line is for the case without Coulomb and surface terms while the dotted line is for a 
larger symmetry energy. The regions of liquid and gas phases and the coexistence are indicated explicitely for the dotted curve. 



loops, the effect of a surface term is to shift the coexistence surface to lower pressure. The loops are still symmetric 
around y = 1/2 and slightly compressed by the droplet 's surface tension. The shifting by the surface term can also be 
seen from other curves with Coulomb terms included. The solid curve contains the symmetry energy, surface energy, 
and Coulomb energy while the dash-dot-dot-dotted curve has symmetry and Coulomb effects. The main effects of the 
Coulomb term are: 1) to shift the intersection point of the two loops to ue = 0.41 for the parameters considered and 
to higher pressure and 2) to make the loops very asymmetric with the coexistence loop for y < he becoming smaller 
and the coexistence loop for y> De becoming larger. Most of the loop for y > yE was not studied and thus was not 
closed in previous studies 0,13 due to the restriction of y < 1/2. The existence of coexistence pair with more proton 
in gas than in liquid was first pointed out in Ref.0 but in a small window of pressure since we restriced y < 1/2 
which allow only up to y — 1/2 for the gas phase. For each case, two loops intersect at yE which occurs at the lowest 
pressure while the maximum pressure for each loop has a different y value. The liquid and gas coexist inside of these 
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two loops since the boundaries are the pure liquid and gas phases. The liquid and gas phases are separated outside of 
these loops; liquid phase for y value between two maxima and gas phase for y value outside the two maxima. These 
are indicated explicitly for the dotted curve. A system with y < ys inside the coexistence region is composed of a 
proton richer liquid and a proton deficit gas compared to the original y while a system with y > yE in the coexistence 
region has more protons in gas phase than in liquid component. The Coulomb term shrinks the coexistence region 
for y < yE while it enlarges the coexistence region for y > yE (Fig^}. 

The boundaries for mechanical and chemical instability are easy to generate. Of these two instabilities, the region 
of mechanical instability is somewhat easier to obtain and follows from (dP/dp) Vt T — which leads to a quadratic 



equation, ap + 2bp 



0. Using Eq. (J2J) , this condition gives 
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where d — t^^-- The presence of d represents a first order quntal correction to the non-degenerate limit. The 
Coulomb term and symmetry energy terms appear in b. The critical point (y c (T) , P C (T)) or (P c (y), T c (y)) occurs 
when b 2 = ac with the critical density p c — —b/a. 

We note that the coefficient b can be rewritten in a form involving b = di(y—yE) 2 +d2 with d\ = (|(| + ^o) a o + d) + 
C and d 2 = (l + C/[(|)(| + xo)&o + d]) — ( a o — f )■ Thus b is symmetric around yE and since y does not appear 
in a or c, the solution p± as a function of y is also symmetric around yE- Therefore, the mechanical instability curve 
of P versus p will have a symmetry that results in one curve instead of a loop. Also the mechanical instability loop in 
P versus y and y versus p are symmetric around yE- This fact is not true for chemical instability and the coexistence 
curve. These facts will be shown in Fig. |21 

The boundary of chemical instability can be obtained from (dp q /dy)p_T — 0. These conditions for proton and 
neutron give the same relation since ydp, p +(l—y)dp n = (1/ p)dP. For two component system with a fixed total fraction, 
if one component changes from one phase to other then the other component must compensate this concentration 
change since they are not independent of each other. The chemical instability condition can be rewritten in terms of 
derivatives of chemical potential and pressure with respect to the variables p and y which appear in Eq. J2J) and Eq. © . 
Namely, the chemical instability boundary can be obtained from (dP/dp) Vi T — {dP/dy) Pj T(dfiq/dp)y t T/(dfJ, q /dy) Pt T- 
Then using the equation of state and chemical potential equation, the chemical instability boundary is determined by 
solutions to a cubic equation 



C 3P + C2P + c\p + c = 



(6) 
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32(- +x )a + 12C+ 12d 
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We note that the term asymmetric in y «-> (1 — y) is the coefficient ci, while ci and C3 are invariant under this 
interchange. The solution to the cubic has three roots. For a physical situation, two roots will be positive and 
physical and the third negative and unphysical. Since c and c 3 are positive, the product of the 3 roots is negative. 
For each y, the solution to the cubic generates a chemical instability curve in p and T and a three dimensional surface 
in p, y, T. Each point on this surface has corresponding values of P and p q . Cutting the (p, y, T) surface of chemical 
instability with a fixed T plane gives a p, y curve for the chemical instability and on this curve P and p q are specified. 
Then using this fact, we can generate p q versus y, P versus p, and P versus y curves for the chemical instability shown 
in Fig[3 Similar mechanical instability and coexistence surfaces can be generated. These results are also shown in 
FigH 

Figure [21 contains plots of y versus p, the loops of p p and p n versus y, curves of P versus p and curves of P versus 
y. Interesting points occur at dp/dy = on these curves which correspond to the cusp points at yE in P versus p 
which have two different branches for the coexistence curve and chemical instability curve. The mechanical instability 
region has only one curve, which starts and ends at yE, because of the symmetry of it around yE- The upper curve 
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FIG. 2: Figures show the coexistence curves (dashed line), chemical instability curves (thick solid line) and mechanical 
instability curves (thin solid line). Also shown are the dp q /dp = curves for proton (dotted line) and for neutron (dash-dotted 
line) at T = 10 MeV. The Coulomb and surface terms are included here. 



corresponds to the proton richer side with respect to y& and the lower branch is for the neutron richer side of yE- 
These two curves form a closed loop in a constant T plane. The upper branch did not appear in Ref.0] due to the 
restriction of y < 1/2. When the Coulomb force is turned off these two branches in the coexistence curve and chemical 
instability curve would collapse into one curve also due to the symmetry around y = 1/2. The maximal asymmetry 
points correspond to dy/dp — 0, one at low y and one at high y. In the y versus p curve, the mechanical instability 
loop is inside the chemical instability loop except at two points where it is tangent to it. The tangent is at y&. The 
chemical instability loop in y versus p is inside the coexistence loop and tangent to it at two points, one at y < yE 
and one at y > yE- These values of y correspond to the critical y c (T), points with maximum pressure of each branch 
having the same densities both in liquid and gas phase, at a fixed temperature, 10 MeV for Fig|5] For P versus y the 
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mechanical instability loop is symmetric about tje, while the chemical and coexistence curves are not when Coulomb 
forces are present, but are symmetric when Coulomb forces are turned off. The Coulomb instability leads to proton 
diffusion above ue and neutron diffusion below yE- The chemical potential p, q and the dp, q /dp = curves show that 
the role between proton and neutron is exchanged around yE- In Ref.Q, the dp, p /dp = curve had two disconnected 
portions due to the cut at y = 1/2 and was quite different from the dp n /dp = curve. 

Similar calculations as shown in FigElfor the cases without Coulomb and/or surface terms can also be done even 
though we did not included the figure in this paper to keep the fiqure simple. The effects of each term when treated 
separately are as follows. The surface energy brings the pressure down, increase the chemical potential of the upper 
part of each loop, and slightly enlarges the loops in p-y plot. The Coulomb term breaks the y <-> (1 — y) symmetry 
and shrink the p-y plot to the lower density and proton concentration while increasing the lower part of the chemical 
potential and pressure curves. 

In summary, in this paper we investigated properties of the coexistence region, and mechanical and chemical 
instability regions when asymmetric Coulomb forces are also included besides volume, surface, and symmetry energy 
terms in a mean field description of hot nuclear matter for the whole range of proton fraction < y < 1 (only the 
range of0<y<l/2 was considered previously for the coexistence region 0, Q and for the instability region 0]). 
Including the Coulomb interactions leads to a changes in many features that are physically important. These include 
a shift in yE, the value of the proton fraction at which the liquid and gas phases would coexist at the same proton 
fraction from y = 1/2 to a value of y closer to the valley of /3-stability. The shifting of yE from 1/2 to a lower 
value by the addition of the Coulomb force is an important effect for isospin fractionization or distillation 9] . An 
investigation of the chemical instability properties shows that proton diffusion will occur above the value of yE while 
neutron diffusion occurs below it. By extending our calculation to values of y > 1/2 > yE we were also able to show 
some new features of the coexistence and instability regions. First, the constant T cuts of the corresponding surfaces 
now form closed loops while the constant T cuts of d(i q /dp — surfaces form continuous curves. Previously, the 
proton curve dfi p /dp = had two disconnected pieces due to the restiction y < 1/2 3]. The main behavior of the 
temperature dependence of these surfaces remain the same as in Ref . 0] with the role of proton and neutron exchanged 
in going from the y < yE region to the y > yE region. The mechanical instability region in P versus p or P versus 
y is symmetric around the shifted yE when Coulomb interactions are included in our model. This symmetry does 
not occur for the coexistence curve and chemical equilibrium curves. Rather, these have new branches which appear 
in P versus p and the two branches cusp at the two p's which have y — yE- The upper branch corresponds to the 
proton rich side of yE, or y > yE, while the lower branch corresponds to the neutron rich side of yE, or y < yE- 
The coexistence curve now has two asymmetric loops in P versus y that meet at the shifted yE- The upper branch 
of the instability curve with y > 1/2 and coexistence loops at y > yE were missing in previous studies. For systems 
with y > yE, the gas phase can become proton richer from proton diffusion induced by a Coulomb instability of the 
liquid phase. For y < yE, the gas phase will be neutron richer from neutron diffusion. This paper also shows how the 
Coulomb interaction breaks some of the symmetries associated with isospin in a hot nuclear system and its associated 
liquid-gas phase transition. The effects of a surface energy terms are discussed. The surface tension associated with 
the presence of a liquid drop shifts the coexistence curve to lower pressure. 
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